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The Keldysh's theory of superfluidity of rarefied electron-hole gas is generalized to a case of 
possible pair polarizability. It was shown that inhomogeneity of the system leads to dipole moment 
which is proportional to the density gradient. The dipole moment appears also near boundaries 
of the system. It was determined that quantized vortices in a magnetic field carry a real electric 
charge. In He II at _H" = 10 T and helium rotation velocity 10^ s~^ the charge density is about 
lO'^ecm"'', where e is the electron charge. 



PACS numbers: 67.90.+Z 

Recent experiments [H - [1] revealed that motion of 
superfluid helium is accompanied by appearance of elec- 
tric fields in the system, and the effect disappears at the 
transition to the normal state. These experiments at- 
tracted a great interest and stimulated a number of the- 
oretical investigations Q - @ in which attempts were 
made to explain these phenomena. However, the physi- 
cal nature of the effects observed is not understood until 
the present. In this situation it seems to us that for a 
qualitative explanation of the electric effects in a mov- 
ing superfluid system it is reasonable to proceed from a 
simple model which allows a consistent microscopical so- 
lution. Our theory may be considered as a purely model 
one. But it can also quantitatively describe two real 
physical systems: spin polarized atomic hydrogen and 
electron- hole gas (exciton gas). Below we will limit our- 
selves to a model of rarefied electron-hole gas in which 
the size of a pair is much less than the average distance 
between the pairs. The properties of this gas were stud- 
ied by Keldysh [lo'|, who showed that the behavior of 
this system in the superfluid state can be described by a 
complex order parameter <&(ri, r2), which has a meaning 
of the pair's wave function. Ignoring the pair's internal 
structure (i. e. assuming ri — r2), Keldysh obtained an 
equation for $ which coincides with a Gross - Pitaevskij 
equation. As far as the effects interesting us can appear 
only if the electron and hole coordinates do not coincide, 
it is necessary to find the function <i>(ri,r2) at ri r2. 
This problem is solved in our work. We obtain an equa- 
tion for <i> at ri ^ v^. It turns out to be a nonlinear 
integro-differential equation which can be solved if the 
order parameter varies slowly at lengths of the order of 
pair's size. Using the obtained $(ri,r2) we can find the 
superfluid density Ps (r) , the superfluid velocity Vg (r) and 
the density of the system's dipole moment P(r). 

Let's proceed to the solution. We will consider 
electron-hole gas with particles which interact accord- 
ing to the Coulomb law. Due to the attraction between 
electrons and holes they form bound pairs, and in the 
spatially homogeneous case the ground state wave func- 
tion coincides with the BCS function. However, we are 



interested in the spatially inhomogeneous case, when the 
density and the velocity of pairs can be functions of co- 
ordinates and time. Keldysh [To'| argued that in this case 
the wave function |<I>) can be looked for in the form 

1$) = expj-y" [*+(ri)$(ri,r2)e^^**+(r2)- 

- h.c. dridrsjlO). (1) 

Here ^'e(r) and ^'h(r) are Fermi operators of destruction 
of electrons and holes correspondingly at a point r. It is 
useful to mention that in the spatially homogeneous case 
Keldysh's wave function transforms to the BCS function, 
in which Wk = cos<I'(k), = sin$(k), where $(k) are 
Fourier components of $(ri, r2). The unknown function 
<I>(ri, r2) can be found varying the system's energy (more 
exactly, the functional $}) by Calculating the 

functional $} in the low density limit, varying the 

difference E — /iiV, where N is the total number of parti- 
cles, and equating the result to zero, we obtain an equa- 
tion ($+(r,r') = $*(r',r)): 
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2mh dY\ 2me dr^ 
+ j i?(ri,r2,r3,r4)$(ri,r3)$+(r3,r4)x 

X $(r4, r2)dr3dr4 = /x$(ri, ra). (2) 
Here the kernel R equals to 
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Constructing an approximate solution, we must take 
into account three characteristic energies in our prob- 
lem: pair binding energy eq, energy of interaction be- 
tween the pairs (per one pair) tint = gn, where g is 
the interaction constant (we will calculate it below), 
n is the density of the pairs, and energy depending 
on inhomogeneity einhom ~ /2ML^, where L is the 
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characteristic scale of the inhomogeneity. We will sup- 
pose that these energies satisfy the following inequali- 
ties: eo ^ Cm* > tinhorn- The first of these inequali- 
ties is a consequence of low density of pairs n supposed 
above. The second inequality means that we consider 
only weakly inhomogeneous states, where the character- 
istic inhomogeneity scale L is greater than the coherence 
length ^ = {fi? /2MgnY/'^ . As far as the energy eq is asso- 
ciated with motion of the electron and the hole relatively 
to each other, and other energies represent motion of the 
pair as a whole, the solution of ^ can be looked for as 
$ = [0(0) (ri2) + 0W (1-12; Ri2)]*(Ri2) {ri2 = ri - rs, 
R12 is the coordinate of the pair's center of mass), and 
0(^) is eint/^o times less than The ^ function (order 
parameter) describes motion of the pair as a whole. 
In the zero order approximation 



(4) 



Here Hq is the Coulomb problem Hamiltonian, /iq = Eq, 
= where ep and </)o are correspondingly the 
ground state energy and wave function. 
In the next approximation 
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2M dUl^ ' 
i?(ri, r2, ra, r4)0o(ri3)0o(r34)0o(r42) x 

X 4'(Ri3)**(R43)*(R42)rfr3dr4. (5) 



Equation ([5]) has a solution only if its inhomogeneous 
part is orthogonal to the solution of the corresponding 
homogeneous equation, i. e. 0o- The solvability condi- 
tion leads to the following equation for ^P: 



Here the interaction constant g equals to 



0. 



(6) 



J R(ri,r2,r3,r4)0o(ri3)0o(r34)'/'o(i"42)0o(ri2)x 

13 

X (ir3dr4dri2 = —ne'^a^. (7) 

While solving the equation ([5]) we must take into account 
that in the last summand, due to existence of three 00 
functions, arguments of all functions must be close to 
R12. Solving then ([5]) we find that the part of ^^^^ which 
makes a nonzero contribution to the dipole moment is 
equal to 



(1)= J dr[,Gir,„r[,)Fir'M, 



g|*0(Rl2)| 

9Ri2 



(8) 



Here G is the Green function of the Coulomb problem, 
i<"(ri2) is an even function of ri2 which explicit expression 
is not shown here because it is too cumbersome. 



Now we can find the dipole moment density of the 
system. 

P = y" [0(°)*(r) + 0(^)*(r)] **(R) X 

X er U(°) (r) + ^^^^ (r)j *(R)dr. (9) 



Substituting here the 0o a-nd 0i functions, we obtain 
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Here 7 = A is a numerical coefhcient, A « 120. 

We see that an inhomogeneity in the system leads to 
appearance of the dipole moment. It becomes zero if 
nih — trie, i. e. the dipole moment is caused by the 
mass asymmetry of the pairs. Usually mh > me- In this 
case in an inhomogeneous system the dipole moment's 
positive end points opposite to the density gradient. 

A real physical system is always limited with vessel 
walls. Interaction between atoms of the system and ves- 
sel walls causes a surface dipole moment. The dielectric 
constant of the electron- hole gas is denoted by ei , of the 
vessel walls £2- The potential of interaction between the 
pair and the surface is found by means of electric image 
method. At large distances to the surface Z ^ z ~ as 
{Z > 0) the potential equals to 
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(11) 



This potential appears in the main equation ([2]) and 
changes both 0(r) and ^'(R). Solving the equation for 0o 
using the variational method, we find the surface dipole 
moment, which is directed normally to the surface. 
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We again see that the dipole moment becomes zero if 
nih — rrie, and when rrih > me (7 > 0) its positive 
end is directed toward the medium with lower dielectric 
constant. Near a metal surface (£2 — > 00) at rrih > rrie 
the dipole moment Pz > 0, i. e. the electron is closer to 
the surface than the hole. Near a boundary with vacuum 
(£2 = 1) the opposite case is realized. 

The potential pTjl influences also the order parameter. 
The equation for it obtains the form 



1 £2 



2£i£2-h£l Z3 



(13) 

Outside the immediate vicinity of the boundary the first 
summand in (|13p can be omitted (Thomas - Fermi ap- 
proximation). Consequently we find that the addition 
to the order parameter caused by Van der Waals interac- 

tion with the surface, has the form I^E'I = t^^^t^^s^- 
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This inhomogeneous addition, according to ([T0|) , leads to 
an extra dipolc moment 



P. 



157— ^^—^en^ 
ei £2 + ei Z'^ 



(14) 



This expression must be added to (|T2l) . 

To understand more clearly the cause of the dipolc mo- 
ment in the multiparticle system considered, it is useful 
to utilize the results for a system consisting only of two 
identical atoms. The Van der Waals interaction between 
the atoms leads to the shift of the electron density distri- 
bution center of each atom towards the opposite atom, 
and dipole moments appear (see e. g. (ll|V 



After multiplying by n this expression coincides in letters 
with (dni). Thus, knowing the electron density distribu- 
tion in a system of two atoms ()15|) . we can qualitatively 
reproduce the results of the microscopical calculation. 

Only stationary states were considered until now. A 
more general state is the Keldysh wave function where 
$(Ri,R2) is replaced by <i>(Ri, R2, t), and its depen- 
dence on t is slow compared to exp(-|/it). As the result, 
right hand sides of equations ^ and © will acquire an 
extra summand ifi^ ■ 

It is interesting to find out the relationship between ob- 
tained results and an affirmation stated by Melnikovsky 
that acceleration of dielectric leads to its polarization, 
and 



Pi = (Pi 



-P2 



-(P2 ■ n) 
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Ri2/J^i2- In particular, for two helium atoms 



18.4ea|. 



Here n 

In a rarefied system the dipole moment of an atom is 
formed by adding the dipole moments which appear in 
this atom as a result of interaction with all other atoms. 
As follows from (fTS]) . in a hypothetical medium where at 
Z > the atomic density equals to ni and at Z < it 
equals to 712, the z component of the dipole moment of 
an atom at the point (0, 0, Zq) {Zq > 0) equals to 



p,iZa) = -D J 



[nie{Z) + n20{-Z)]{Z-Zo) 
[p2 + (Z - ZoYV 



(fpdZ. 



(16) 

Here 9{Z) is the Heaviside step function. When writing 
this expression we took into account that the z compo- 
nent of the dipole moment of a pair of atoms is obtained 
multiplying psp by cosine of the angle between the z axis 
and the vector connecting these atoms. As far as the dis- 
tance between the atoms cannot be less than the size of 
an atom, p and Z must be cut off at the lower limit at 
the Bohr radius as- Further we assume Zq ^ as- After 
a simple integration in (I16p we obtain 



e-1 1 , 5v 
M— . 

47r 2Ze dt 



(19) 



Here (and only here) Z is the atomic number. Below 
Z — 1. Denoting as C/(R) the interaction energy of the 
atom in the point R with all other atoms and taking into 
account again that for a rarefied medium e = 1 + Airna, 
we obtain from 



P(R) 



Q;n(R) 
2^ 



dU 



(20) 



If (not very justified) we spread this expression to two 
atoms at points R and R', then at |R — R'| ^ qb the 
Van der Waals attraction forces will act between these 
atoms and U = ~C/\R - R'|'', where C « e^a|. As 
far as P = np, (|20p yields an expression for the dipole 
moments of a pair of atoms which coincides by sign and 
by order of magnitude with psp . This result confirms 
applicability of (|19p down to atomic length scales. 

In the case of a dielectric medium the total dipole mo- 
ment of the system equals to 



PdV 



a 
2^ 



(9 77 

nUdS - I U—dV 



(21) 



TT 712 - "1 

Pz = -zD- 



(17) 



If we take into consideration that for a rarefied system 
ei_2 = 1 + 47rani_2 and the polarizability a — Ca^, 
where C is a numerical coefficient (for hydrogen C = 
for helium C = j), it is easy to prove that coincides 
in letters with ((17)) multiplied by n to obtain the dipole 
moment density P. 

Above we assumed the densities rti^2 to be homoge- 
neous. But near the vessel walls due to the Van der 
Waals interaction the density acquires an addition ~ , 
where Z is the distance from the wall. Supposing that ni 
is a slow function of Z and writing it as ni — uiq + ^Z, 
after the integration in (jl6p we obtain the addition to pz 



^Pz 



TT dn 

64 az' 



(18) 



We have used the expression (PO)) and integrated its right 
hand side by parts. As it was shown above, the pair 
interaction energy (per one pair) equals to gn. Using ([2T|) 
and taking into account U = gn we obtain the volume 
part of the dipole moment 
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P = — —naea\n-^^. 



(22) 



This expression coincides in letters with (jlOp when nih ^ 
TOg. The last inequality is implicitly supposed to be ful- 
filled in 

In the general case the density can be a complicated 
function of coordinates and time. But a superfluid sys- 
tem can possess "characteristic configurations" which are 
of special interest, in particular, rectilinear vortices and 
vortex rings. In the case of a rectilinear vortex the den- 
sity becomes a function of distance p to the vortex axis. 
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As the result a radially directed dipole moment appears 
in the system. 



3A nag p 
IStt e Mp'^ p 



(23) 



This result (with an opposite sign) was first obtained by 
Natsik [3] phenomenologically. 

A more complicated problem about the dipole mo- 
ment of a vortex ring can be solved if we take into ac- 
count that far from the vortex core a relation is fulfilled: 
V(Mw^/2) = —gVn. As far as the velocity field of a vor- 
tex ring is known , it is easy to find the dipole moment 
of the ring as a function of distance to its axis. We can 
also find the total dipole moment of the ring and show 
that it is zero unlike the affirmations found in literature 
that it is nonzero. 

Until now we assumed that there is no external mag- 
netic field. Presence of a magnetic field leads to a cardinal 
change in the situation. The effect is caused by appear- 
ance of a summand f (v x H) • ri2 in the Hamiltonian. 
Here v = is the velocity of the pairs. The system's 

response to this addition is similar to the response to an 
external electric field E^f f — ^{v x H) . Like a real elec- 
tric field £'e/ / leads to appearance of a dipole moment in 
the system: 



(v X H)ri, 



(24) 



where a is the pair polarizability. This dipole moment is 
connected with the Lorentz force acting on the positive 
and negative charges in opposite directions. If we substi- 
tute the velocity field of a rectilinear vortex v in ([M]) we 
will see that a radial dipole moment, as without a mag- 
netic field, will appear around the vortex axis. But it is 
extremely important that this dipole moment, similarly 
to the velocity field, decreases with distance from the 
axis as 1/p. As the result the total surface polarization 
charge J P ■ dS does not depend on p and on the surface 
shape if the surface remains cylindrical. Due to electric 
neutrality of the system this result leads that the vortex 
core possesses the opposite sign charge which equals per 
unit length (H being parallel to the cylinder axis) 



ah 

± 2nHn. 

Mc 



(25) 



The charge sign depends on the sign of vortex circulation. 
If the vessel with helium rotates with frequency Q, the 
density of the vortices equals to = Mil/iTh, and the 
total charge of the vortices per unit volume is 



2an „ 

± — Hn. 

c 



(26) 



For He II in a magnetic field 10 T and at rotation speed 
10^ s~^ the charge density is approximately lO^ecm"'^, 
where e is the electron charge. 

Appearance of electric charge in the vortex core will 
lead to nonzero dipole moment of vortex-antivortex pairs 
(in a 2D system) and vortex rings (in a 3D system). 



We must mention that electric polarization and electric 
charge in the vortex cores in ^He were considered in [l3{ . 
Unlike our work, the polarization and the charge of the 
vortices in '^He are caused by fiexoelectric effect. 

Now let us show that due to existence of the dipole mo- 
ment P near the metal boundary a second sound wave in- 
cident to it will lead to an oscillating potential difference 
between the metal and the helium. Since the induction 
D = E + AttP must be continuous and in metal D = 0, 
in helium filling the metal vessel an electric field appears 
in the vicinity of the vessel walls: E — —AttP. This field, 
similar to the dipole moment P, is constant at constant 
density n. But when a second sound wave propagates,the 
density n and the dielectric constant e oscillate in time, 
more exactly, they obtain oscillating additions (see be- 
low). As far as D = eE and SB = (5eE + eSE = 0, 
the oscillating addition de leads to similar addition to 
the field S'E — — — E. This addition is the cause of the 
potential difference between metal and helium. 

The oscillating addition Se can be found taking into ac- 
count that the Clausius - Mossotti relation = Aima 
leads to 5e{t) = \{e-\){e + 2)^§^ST{t). 

The expansion coefficient of liquid helium /3 = 
can be taken from the experiments (see [3|). For ex- 
ample, /3 = -0.1 • 10-3 K-i at r = 1.138 K and 
(3 = -26.34 • 10-^ K-^ at T = 2.15 K. Since e - 1 = 
5.7-10-^, at T = 2.15 K the oscillating addition equals to 
6e — —1.5 • IQ-^^T where ST is in kelvins. The constant 
field at the surface of helium substantially depends on 
zq. Choosing zq — 5 ■ 10* cm, we obtain E ~ lO^V/cm. 
In this case the oscillating potential difference between 
metal and helium (50 = SeEzo at T = 2.15 K will be 
equal to S(j) « —5 • IQ-^JTV. This value is close to the 
experimental one, but strongly depends on temperature. 
At T = 1.138 K we obtain Scf) « lO-^f^TV, i. e. two 
orders less than experimental data. 

Comparing theory and experiments we must take into 
consideration that without special preparing of the ves- 
sel containing helium its walls are covered with several 
layers of adsorbed atoms. The interaction of these atoms 
with metal will induce dipole moments in them and elec- 
tric fields. The polarizability of adsorbed atoms is almost 
two orders greater than the helium polarizability, there- 
fore, the induced electric field in the adsorbed atoms will 
be correspondingly greater. Estimates given above allow 
to expect that the scale of the effect will be the same as 
in experiments. As far as the superfiuid transition tem- 
perature is not an exceptional point for adsorbed atoms, 
there are no reasons to expect a substantial temperature 
dependence of S(j) in the temperature interval 1.5 — 2K 
in which experiments were performed. 

It should be noted that the effects predicted in this 
work do not give a comprehensive explanation of the ex- 
periments In our opinion, for such an explanation 
we need a new key idea. 
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